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ABSTRACT 

The Scanning Multichannel Microwave Radiometer (SMMR) is a 5-frequency ( 6 . 6 . 10.7, 18, 21, 
and 37 011/ dual polari/cd microwave radiometer which was launched in two separate satellites, 
Nimbus-7 and Scasat in l l >78. A formalism is developed which can be used to interpret the data in 
terms of sea surface temperature, sea surface wind speed, and the atmospheric overburden of water 
vapor and liquid water It is shown with reasonable instrumental performance assumptions, these 
parameters can be derived to useful accuracies. Although the algorithms are not derived for use in 
rain; it is shown that, at least, token rain rates can be tolerated without invalidating the retrieved 
geophysical parameters. 
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AN ALGORITHM l()K Rl I Kll VA1 Of (HI AN SURI AC E AND 


A IWOSIMII RIC PAR AM I II RS I ROM till OBSERVATIONS OP Till 
SCANNING MULTIC IIANNI l MICROS AVI RADIOMI I I K (SM MR) 

SECTION I 

INTRODUCTION 

The Scanning Multichannel Microwave Radiometer (SMMR) was carried aboard both the 
Nimbus- 7 and Seasat satellites, both of which were launched in 1478. It is a 5 frequency, dual 
polarized microwave radiometer which scans an 80 cm reflector such that the antenna beam sweeps 
a conical arc of 50 ’ with a cone angle of 42° measured at the spacecraft. This provides constant 
incidence angle of approximately 44 (for an 8(H) km orbit) at the Earth's surface when the cone 
axis is aligned with the local vertical (its nominal orientation). The primary purpose of this instru- 
ment is to measure the sea surface temperature, sea surface wind speed, and the atmospheric content 
of water vapor and liquid water (clouds) under nearly all weather conditions. It is the purpose of 
this paper to derive a formalism for extracting these parameters from the measurements provided by 
the SMMR. A similar derivation was provided by Chang and Wilheit ( 1479) for interpretation of a 
three channel system. We will refer to this paper for computational details and only repeat enough 
for conceptual self-sufficiency in this paper. We will also draw heavily on review articles by Wilheit 
(1978) and Wilheit et al. (1478) for background discussions. 

In Section II. we will treat the relationships between the geophysical variables and the bright- 
ness temperatures. An approach for extracting the geophysical parameters from the brightness 
temperature will be developed in Section III and numbers will be calculated based on the models of 
Section II. In Section IV we will develop an algorithm for approximating rain intensity from the 
radiometer measurements. In Section V, we will examine the sensitivity ol these algorithms to rain 
filling all or part of the field of view . In Appendix I we will give computer programs for all the 
retrievals and explanations of them. 
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SECTION II 


BACKGROUND PHYSICS 

Because of the validity of the Rayleigh Jeans approximation in the microwave region, the law of 
radiative transfer may be expressed quite simply 

^*--T(T,-t„> 0) 

where T B is the microwave radiance expressed as an equivalent black body temperature. 7 is the ab- 
sorption coefficient at the point x and t e is the thermodynamic temperature of the medium at that 
point. 

The derivative is along the direction of propagation. It must oe emphasized that 7 is only the 
attenuation due to loss mechanisms; scattering and reflection must be accounted for by appropriate 
re-distribution of radiation among directions of propagation. 

A typical Earth viewing geometry is shown in Figure I . In order to calculate the brightness 
temperature expected for such a situation one first divides the atmosphere up into a number of 
layers each of which may be characterized by a uniform temperature and absorption coefficient. 
Beginning with the 3°K. cosmic background, the radiation is then propagated from the top down- 
wards through each layer according to Equation I . A fraction R (the reflectivity) is reflected from 
the surface and the surface also radiates an amount Ft where E is a quantity characteristic of the sur- 
face. the emissivity, having a value between 0 and I and t is the thermodynamic temperature of the 
surface. Arguments based on thermodynamic equilibrium show that F and R are related by: 

F + R = 1 (2) 

Since water is a polar molecule it has a very large dielectric constant at microwave frequencies. 
This results in a large reflectivity (low emissivity) for a liquid water surface such as the ocean. In 
Figure 2. the emissivity. calculated for a smooth water surface using the dielectric constant data for 
20° (' of Lane and Saxton (I952)and the Fresnel relations (Jackson l%2), is shown. Note that the 
emissivity increases somewhat with frequency over the range, and that for viewing at an oblique 
angle, the horizontal polarization yields a lower and vertical a higher emissivity than that for nadir 


viewing where the distinction between the two polarizations breaks down. The emissivity was cal- 
culated both for fresh water and for a 3.5 ‘a NaCJ solution which is a reasonable approximation to 
sea water: there is no consequential difference between fresh and sea water except at the very low 
frequency end of the spectrum, below about 5 C»Mz. Since most solid surfaces have emissivities in 
the range 0.8 to 0.95 and 0.5 is a typical value for water, there is a striking contrast between liquid 
water surfaces such as lakes, rivers, and the oceans and solid surfaces such as land and sea ice. This 
low emissivity of the sea surface also provides a good background for viewing the intervening 
atmosphere. 

The emissivities shown in Figure 2 were calculated for a smooth water surface. However, when 
the wind blows across the surface of the ocean, it generates roughness and foam. Some care is re- 
quired in the quantitative definition of wind speed at the surface. There is a substantial vertical 
gradient in the wind speed in the boundary layer. A turbulent wind shear model for the marine 
boundary layer has been developed by Cardone (1969). In this model, the details of the w ind shear 
depend on the wind speed and on the air/sea temperature difference. In general, there is more shear 
if the ocean is colder than the air (stable) than if the ocean is warmer than the air (unstable). Wher- 
ever possible wind speeds used in the derivation of the emissivity models defined as follows: The 
wind speed is measured at any altitude within the marine boundary layer. Using the actual air and 
sea temperatures, the friction velocity (generally termed *‘U-Star" (U*)) is calculated with the Car- 
done ( 1969) model. Then, assuming neutral stability (air/sea temperature equal), the wind speed at 
20 m altitude is calculated I his last step is merely a translation of U* into a more familiar quantity. 
The results of Wilhcit (i978 II) suggest that this is an appropriate definition. The effect of roughness 
on the emissivity of the water surface depends strongly on the viewing angle and polarization. For 
viewing directly at nadir where the distinction between vertical and horizontal polarizations dis- 
appears, roughness has very little effect on emissivity; for viewing in horizontal polarization at angles 
away from nadir, the roughness increases the emissivity and for vertical polarization for view angles 
between 0° (nadir) and 55° the roughness increases the emissivity while it decreases the emissivity 
for view angles greater than about 55". <s was determined experimentally ny Huliiugu ( 197 i ). Foam 
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cover, on the other hand, always increases the emissivity of the surface. Figure 3 shows the increase 
in brightness temperature caused by a given wind speed as found by Nordberg et at. ( 1971 ) for a 
frequency of 1 9.35 GH. and for nadir viewing. 1 hey found no effect, in their case, for wind speeds 
less than 7 in s and an increase of around I Ks'm above that. This relationship has been investigated 
further b> Webster et al. 1 1975). They examined a frequency range of from 1 .4 to 37 GHz and a 
\iew angle of 38° (both > olurizations). Figure 4 summarizes their results. Here the rate of increase 
of brightness temperature with increasing wind speed is shown for all frequencies observed. Note 
first that viewing in horizontal polarization enhances and viewing in vertical polarization decreases 
the strength of the effect vis-a-vis nadir viewing. Tins is because in horizontal polarization the effect 
of roughness and of foam both increase with increasing wind speed and because viewing in horizontal 
polarization decreases the emissivity of the undisturbed surface thereby increasing the emissivity 
contrast between the foam and open water with respect to nadir viewing. ITie converse argument 
for foam and the relative weakness of the roughness effect account for the vertical polarization case. 
Note also that the spectrum of the net effect is only weakly frequency dependent above about 
10 GHz but it decreases somewhat at lower frequencies. 

W illicit ( I97X II) has examined the roughness effect isolated from the atmosphere and foam 
effects both in satellite data (Nimbus-6 hlectrically Scanned Microwave Radiometer) and in the data 
of Hollingcr (1971). 

It was found that the surface could be treated as an ensemble of flat surfaces and the emissivity 
calculated for each member of the ensemble as though it were an infinite plane surface. The rough 
surface emissivity is then the ensemble average of the individual emissivities. It was found that the 
distribution of surface slopes given by Cox and Munk ( 1954) overestimated the roughness effect 
below 35 GHz, but was acceptable at higher frequencies. In particular it was found ihat the variance 
of surface slopes o’ required to account for the microwave observations is given by: 

o 2 * (0.3 +0.02 f (GHz)) f£35C.Hz 

o’ = o 2 m f^ 35 GHz (3) 
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where I' in the observation frequency and u* m in the surface slope variance given by Cox and Munk 
(1955). A composite surface model which includes the effect of both roughness and foam was pro- 
posed by Wilheit ( 1979). In this model, the rough surface described is partially obscured by foam 
(for wind speeds greater than 7 in/si such that R. the reflectivity of the surface, is reduced by: 

R n , (?.PI s (I - K) Rroighnus 

only 

k * a( I - 1 1 ,0 ) (w - 7 in/s) w>7m/s 

k * 0 w ^7 m/s (4i 

where f is the observation frequency, w is the wind speed, 0 the view angle, I’ the polarization and 

a = 0.006 S/m 

f D = 7.5 (1Hz (5) 

The constituents of the atmosphere which make consequential contributions to the brightness 
temperatures observed in this spectral region are molecular oxygen water vapor and water droplets. 
Molecular oxygen has a series of resonances in the 50-70 (ill/ region and a single isolated line at 
I l l » (illz. The microwave properties ol molecular oxygen in an atmospheric context have been dis- 
cussed by Meeks and Lilly ( 1963) and by Carter et al. ( 1968). A very' elegant model which accounts 
for the overlapping of the pressure broadened lines has been published by Roscnkrantz ( 1975). 

Since the mixing ratio ot oxygen is substantially constant and the absorption coefficient is only very 
weakly temperature dependent, it contributes only a constant offset to surface observations is the 
140(111/ region. It must be included in calculations using any of the above models but is simple 
enough that it will not be dealt with further here. 

Water vapor has a fairly weak resonance at 22.2 (ill/ and very strong ones at 1 83 (ill/ and 
higher frequencies, the wings of which give significant contributions to the atmospheric opacity at 
all frequencies but particularly above about 10(111/. flic absorption coefficient for water vapoi 
based on Staelin ( 1966) is given in Figure 5 for the 1962 U.S. standard atmosphere (AFCRL, 1965) 
for sea level and 8 km altitude. The absorption vanes as a function of altitude because of pressure 
broadening. However, the bulk of water vapor is found in the bottom few kilometers of the atmo- 
sphere. over which there is little change in the pressure broadening. Thus, variations in the vertical 
distribution of water vapor cause little difficulty in retrieving surface parameter. 
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Hu absorption due to non-raining liquid water clouds is well described by the Rayleigh approx- 
imation (Gunn and I ast 1954) u.ing the dielectric data of 1 ane and Saxton (1952). In Figure 6, we 
show the absorption coefficient for a cloud with a liquid water density of I gm/m 1 at three different 
temperatures (0 t20°C). The absorption coefficient is almost precisely quadratic in frequency 

and varies by about a factor of three with temperature over the range shown. Since the spectrum is 
preserved with this temperature variation, the spectral characteristic may be used to correct for 
clouds even though the absolute amount of liquid watei would be uncertain by a temperature de- 
pendent multiplicative factor. Additional discussion of the microwave properties of molecular 
oxygen, water vapor and clouds along with numerical expressions for calculating absorption coeffi- 
cients can be found in Chang and U illicit ( 1978). 

Rain, on the other hand, provides a much more complicated situation. The particles are large 
enough compared to a wavelength that resonant effects occur and the absorption is no longer simply 
proportional to the net water content. Further, the spectrum of the absorption coefficient depends 
on the drop size distribution and scattering must be considered as well as absorption. A model for 
calculating the upwelling brightnesc temperature from a raining situation is discussed by Wilheit et al. 
(I‘>77). At present, there is no adequate technique for correction of observations in the presence of 
rain. In Section IV we will develop an algorithm lor determining rainfall intensity and in Section V 
we will approximate the rain intensity through which the retrieval of surface parameters may be 
performed within reasonable error limits. 



SI ( I ION III 


Rl IK 1 1 VAL OF GEOPHYSICAL PA RAM I II RS 

Microwave instrument are typically limited in their spatial resolution by diffraction in the an- 
tenna aperture Since in SMMR all frequencies share a common aperture, the spatial resolution with 
which it views the I arth’s surface at a pis en frequency is proportional to the wavelength. The SMMR 
resolution varies from approximately 150 km at 6.6 GHz to 25 km at 37 GHz. If the various frequen- 
cies are to be used together to arrive at surface parameters, they must be reduced to a common basis 
so that all the measurements apply to the same area. One is. therefore, limited in the retrieval of any 
given parameter to the spatial resolution of the lowest frequency used. Since resolution is always to 
be desired, it is necessary to reduce the SMMR data to several bases at various resolutions. The 
present scheme which is indicated in Figure 7 is to reduce SMMR data to four gilds The first has a 
resolution of 150 km and has all five frequencies; it is used for surface temperature determination. 
The grid 2 has a resolution of 90 km and has all the frequencies except 6.6; it is used for wind speed 
retrievals. In principle, it is possible to use the sea surface temperature derived at grid I to arrive at 
an approximate sea surface temperature to use as an additional input parameter in grid 2 for refine- 
ment of the wind speed determination. The retrieval algorithm presently discussed does not make 
use of this but the stiucture exists in both the Nimbus and Seasat data processing schemes to do so; 
it remains a possible algorithm refinement to be used as confidence in the instrument and in the 
simpler retrievals grows. Grid 3 has a resolution of approximately 60 km at the three highest fre- 
quencies and is used for the atmospheric water retrievals. Both the wind speed and sea surface tem- 
perature can be used as additional input parameters in future algorithm versions. The fourth grid has 
only 37GH/ data, but at a resolution of 30 km. It is only used, in ocean applications, for providing 
strucuual information to the rain rate retrieval which depends primarily on the IS Gil/ data. A 
paper by Njoku ( 1979) describes the very difficult process of reducing the SMMR outputs to these 
common bases. 

As can be seen from the previous section, the brightness temperatures which would be observed 
upwclling from the atmosphere depend on many geophysical variables. Moreover, several of the 
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variable* (temperature, humidity, and liquid water content) arc, in themselves, continuous functions 
of height. Therefore, this problem has an infinite number of degrees of freedom. In the most gen- 
eral sense, a retrieval of geophy sical vanables from a finite set of observations could never be per- 
formed We must resort to approximations and constrained solutions if we are to derive useful 
information. 

One technique for retrieving geophysical parameters is the statistical technique which was ap- 
plied. quite successfully, by Waters et al. ( l l >75) to the problem of approximating the atmospheric 
temperature profile from microwave brightness temperature measurements in the moleculai oxygen 
band (X * 5 mm). In this technique, which is schematically represented in Figure 8, they begin with 
a statistical data base consisting of a large set of actual temperature soundings. For each of these 
soundings, they calculate the expected upwelling brightness temperature at all frequencies at which 
measurements are made. The temperatures at various heights are then regressed against the bright- 
ness temperature or functions of the brightness temperature to determine the most probable (if the 
statistics are Gaussian) temperature from each level. 

If there is little information about a given level in the brightness temperatures, then the retrieved 
temperature will tend toward the a priori average of the statistical data set and the residual of the 
regression will not be much smaller than the variance of the a priori data set. It is to be noted that 
information about a given level can occur in a given measurement either directly through the equa- 
tion of radiative transfer or indirectly through correlations inherent to the atmosphere itself. 

This technique may be faulted in that any bias of the a priori data set is introduced into the 
retrieval. V/hcn the brightness temperatures are strongly related in a direct sense to the geophysical 
parameters being retrieved, this would be of little consequence, but where the coupling is through 
statistical correlations among the geophysical observables or when the retrieval offers only marginal 
improvement over the a priori statistics and the retrieved values will tend strongly toward the com- 
monplace situations. When tested on actual data, this may well give average statistics which suggest 
a very good retrieval technique yet give poor results in »he extraordinary situation which is often of 
the most interest We. therefore, have modified this technique somewhat. We generate an artificial 
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data set subtending the approximate expected range in all the geophysical variables at issue and 
force any correlations among them to be zero or at least much smaller than expected in nature. As 
a by-product, this relieves us the necessity of providing a true statistical data base for combinations 
of parameters which have been measured sparsely or not at all. For each member of this ensemble, 
we calculate the ten brightness temperatures (5 frequencies, 2 polarizations) as well as the geophys- 
ical parameters of interest here, sea surface temperature, sea surface wind speed, the column densities 
of water vapor and cloud liquid water and a correction for refraction by atmospheric water vapor 
calculated accord .ig to Bean and Dutton (1966) which is used in interpreting the data from the 
Seasat altimeter 

Although the retrieval problem contains, in principle, an infinite number of degrees of freedom, 
there are four principal degrees of freedom here, sea surface temperature, sea surface wind speed and 
the net atmospheric content of cloud liquid water and water v. nor. One would not expect details 
of the temperature profile of the atmosphere nor the vertical distribution of the two phases of atmo- 
spheric water to have a large effect. In fact, the specific frequency choice for the 21 GHz channel 
being a little off the water vapor resonance at 22.235 GHz was made, based on unpublished calcula- 
tions by the authors, in order to reduce the minor dependence on the vertical distribution of water 
vapor to negligible proportions. Similar calculations by Westwater ( 1978) arrived at a similar 
conclusion. 

The quasi-statistieal data base used here consists of 90 surface models ( 10 wind speeds and 9 
sea surface temperatures - see Table I) and 81 atmospheric models made by combining the 9 hand- 
book atmospheric temperature profiles in Table II (AFCRL 1965) with the 9 cloud models listed in 
Table III. These cloud i '>• • Jels were simply invented for the present purpose. Observations of cloud 
liquid water content arc almost non-exisienl. These were chosen to give a w ide range in the net 
liquid water content and vertical distribution of water. The observations of Wilheit and Fowler 
( 1977) would suggest that these models are. at least, reasonable. T he combination of all 90 surface 
models with all 81 atmospheric models provides a net data set with 7290 composite models. In 
addition tothe brightness temperatures and geophysical parameters another variable must be 



considered Real spacecraft do no always lly in their nominal attitude but may depart from it by 
several degrees hi pitch. roll. and or yaw I Ini- the angle at which the sensor views the Earth's 
surface can vary as a function ol tune and scan position. This effect is included in the calculations 
by alternately using -4H an *. 50" as t lie view angle as the computation steps through the HI different 
atmospheres. I'he view angle is then treated as an additional observable and is included in the re 
giesston like an additional brightness temperature. 

Iliese models subtend the approximate range expected for each of the parameters but the ex- 
pected correlation* are left out. I veil so absurd a combination as an arctic winter atmosphere over 
an ocean with a tempera are is 2‘ ,l >"k is included This would tend to make simulated retvieval re- 
sults somewhat pessimistic. A weak correlation is introduced between cloud liquid water and water 
vapor content in two ways I list, wherever a cloud introduced, the relative humidity is increased to 
100 (with respect to the liquid phase). I his introduces some variability in the vertical distribution 
ol watei vapor ovei that given by the handbook atmospheres which are generally rather smooth. 

I lie second effect included is spontaneous nuclcation t Mason l l >7l V 1 iquid water droplets in 
clouds cannot be super-cooled lower than -40 (' without freezing. Since ice clouds are transparent 
at 1,'tesc frequencies, the cloud liquid water content is set to zero for any level at w hich the tempera- 
tun is less ili.in 2.0 k I '.us reduces the liquid water content for the colder, and thus drier, 
atmospheres. 

I'he regression technique is inherently linear but the problem at hand has some mild non- 
linearities I lie present algorithm uses two techniques for mitigating the e lect of those non- 
linearities I lie Inst has already been alluded to the use of functions of brightness temperature 
rather than the brightness temperatures themselves. One is free to choose any function of brightness 
temperature for use in the regression, the appropriateness of any particular choice is determined by 
the residuals ol the regression analyses On the basis of heuristic arguments by Wilheit et al. ( I ‘>77) 
we have chosen to use the function: 

F(T„) h In (280°K -T„) 

m lieu of brightness temperature for the channels most affected by the atmospheric constituents. IH, 
21. and .47 011/. 
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rho second approach for dealing with non-linearities is to iterate the solution. The application 
of iterative solutions in the present context is illustrated in f igure 9. A set of retrieval matrices are 
made using the entire ensemble and also various restricted : ubsets of the ensemble ol geophysical 
models. I he data are first interpreted using the retrieval based on the entire ensemble and the ap- 
proximate values for the geophysical parameters are then used to select the matrix derived from the 
most appropriate restricted ensemble. In principle, one could procede to even more and more re- 
stricted ensembles. In the present algorithm, the application of this iteration is limited. The non- 
linearity to be addressed is the onset of foam generation at wind speeds of 7 m/s. This produces an 
abrupt change in Ihe slope of brightness temperature ven.is wind speed. The solution for the sea 
surface parameters (temperature and wind speed) are significantly affected by this abrupt change. 
Thus in both cases the solution is iterated once to decide whether the wind speed is . hove or below 
7 m/s. The solution for atmospheric parameters is not particularly improved by such an iteration 
and so the atmospheric retrievals are not iterated. It was also found that the use of the 37 C»H/ 
channel only aided significantly in the case of cloud liquid water. So. to aid in the isolation of pos- 
sible problems with real data, it was not used in any of the other retrievals. 

Ihe assumed noise values were based on performance measurements of the engineering model 
of the SMMR instrument and approximations for the effect of averaging over the various grids. As 
a tatter of judgement, no noise values less Ilian 0.2' k were assumed. Ihe values used are 0.2' k 
for the grid I and 2 retrievals and 0.5 k lot the grid 3 retrievals The assumed uncertainty in the 
incidence angle is 0.2 of angle The resulting equations for the retrieval are as follows. 

Wind Speed Retrieval - V' id Speed Unknown 
WS (m/s) = -4(>5.3 + d>2lo T B 10.7 y 
+ .28 3 I mo 7 H + l(>8. - in (280. - T B) H v ) 

-86.31 In (280. - I B , „ „ ) 

- 15.84 In (280. - T B2 , v ) 

-37.18 In (280. -T B21H ) 

+ 2.357 0, Nt - : (1.8 m/s / 1 2 m/s) (6) 

1 1 



Wind Speed <7 m/s 
WS(m/s) 3 - 523.9 - .2229 T, 


BIO 7V 

+ (’056 T„ l0 7M + 130.3 In (280. - T B18V ) 
-39.19 In (280. - T B , 8H ) 

+ 10.24 In (280. - T B2IV ) 

-32.75 In (280. -T B2JH ) 

+2.999 0 |NC ;( 1.6 m/s / 1 .9 m/s) 


W nul Speed > 7 m/s 
WS (m/s) = -338.4 + .3115 T B)0 7V 
+ .4509 T„ I0 7I( + 151.8 In (280- T„,, v ) 
-9|.!7ln(280.-T B1BI! )-26.66ln(280.-T B21v ) 

I 2.89 In (280. - I 1 ii ) + I 432 fl| NC . ( I .Om/s / lO.Om/s) 

Sea Surface Temperature ( WS < 7 m/s) 

•SST (°KI = -149.1 + 1 .677 T B6AV ♦ 1.666 T (# t „ 

“ •2767 r B | 07V -.5590 r BI07H 

+ 46. 1 7 In ( 280 - T B , 8V ) + 3.097 In ( 280. - T„ , 8H ) 
.9162 In (280. - r B2lv )- 12.54 In (280. -T B21H ) 

- .5850 0| NC (0.7 o C/8.6°C) 


Sea Surface Temperature (WS >7 m/s) 

SST (°K) = 188.9 + 3.040 T B6 6V -I.I88T B6AH 
-.709 I H)07V + .2405 T bio . 7H 
- 6.1 14 In (280. - T BI8V ) + 20.37 In (280. - r BI8H ) 
-4.003 In (280. - T B2IV ) + .986 In (280. -T B21H ) 
-4.735 0, Nr ;(1.5 o C78.6°C) 



Cloud Liquid Water 

CLW <mg/cm 2 ) « 246.1 - 51.72 In (280. - T B1SV ) 

+ 1.14.4 In (280.-T B , 8H ) +46.14 In (280. -T B2IV ) 

+ 24.95 In (280. - r B21H )- 155.5 In (280. - T„ J7V ) 

-36.63 In (280. - T BJ7H ) - 3.391 0 tNC ; (4.5 ing/cm 2 / 23 mg/cm 2 ) (II) 

Water Vapor 

WV (gm/cm 2 ) « - 9.784 

+ 6.927 In (280. - T B1HV ) + 5.361 In (280. - T B1(iH ) 

-4.518 In (280.-T B21v ) -6.081 In (280. -T B21M ) 

+ .039 0, NC (.15 gm/cm2 / 1.7 gin/cm 2 ) (12) 

Altimeter Wet Correction 
PL (cm) = -51.95 + 34.37 In (280. -T B18V ) 

+ 37.1 5 In (280. - T B18H ) - 22.64 In (280. - T B21V ) 

- 39.70 In (280. -T b ,,„) - .18 0 , nc ;(.9cm/10cm) (13) 

Hie numbers in parenthesis following each equation represent respectively the RMS residual 
error of the regression and the RMS variability of the parameter in the original data set. Fortran IV 
programs which execute these equations and iterations are given in Appendix I. 
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si chon iv 


RAIN AKIORITIIM 

A model lor quanliialivcly relating microwave brightness temperature to rain rati over ocean 
area was developed by Wilheit et al. ( I ‘>77). In their model, the rain is treated as a Marshall-Palmer 
( 1‘>4X) distribution of water droplets extending from the surface to the freezing level (0 l T isotherm) 
with appropriate water vapor and non-raining cloud distributions. It is assumed that the water above 
the freezing level is frozen and thus essentially invisible to the frequencies of interest here. The full 
Mie ( l‘>OXt scattering and absorption properties of the ensemble of droplets are calculated and the 
equation of radiative transfer solved with iteration to account for scatter.. '«?. The problem is domi- 
nated by the microwave absorption by the rain w ith the scattering providing small corrections. 

There are two major problems in the application of this model to the i. terpretation ol measure- 
ments The first is spatial resolution. The relationship between brightness temperature and rain rate 
which results from this model is highly non-linear Thus if the rain intensity varies over the radiom- 
eter’s instantaneous Held of view, the brightness temperature which the antenna averages over the 
field of view in a linear manner will not yield an average rain rate but. in general, a lesser rain rate. 

V the 25 km resolution of the Nimbus-5 Electrically Scanned Microwave Radiometer (ESMR). this 
problem is serious; at the (»0 km resolution of the IX f ill/ channels of SMMR. the problem will be 
much worse. The 57 (ill/ channels have a resolution of 50 km but at rain rates on the order of 
4 mm lu or greater, the atmosphere is too opaque to yield useful rain intensity estimates. This rep- 
resents too low a dynamic range to be useful. The I8(illz can measure rain rates up to about 20 
mm hr. a much more reasonable dynamic range. 

The second problem is that knowledge of the freezing level is required. Raoet al. (l‘>77), m 
producing an atlas of oceanic rainfall from the Nimbus-5 TSMR data, used climatology to provide to 
the nearest I km a freezing level based on location and season. Since the freezing level was given in 
I km increments, changes in the freezing level produced abrupt changes in rain rates. A more serious 
problem was found at high altitudes where the climatological freezing levels were low (1-2 km). It 
was found that these low freezing levels caused the rainfall determinations to be clearly excessive. 
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As an ad hoc solution. I he free/ing level was set to 3 km for these cases. Since the maximum value 


for freezing level is about 5 km it was decided th.it the error contribution would not be too severe if 
the assumed freezing level were fixed at 4 km for the present. I Ins can cause errors of ca. 25'/ in re- 
trieved rain rate; the resolution problem can cause larger errors. The primary purpose of this rain 
algorithm is to Hag the retrievals when the rain rate is excessive rather than to provide ram data 
per se. The effect of the rain on any channel depends first on the net opacity of the rain column at 
that frequency so fixing the freezing level should not have too great an effect on rain rates derived 
for this purpose. A more sophisticated rain retrieval should be used to derive rain data lor climato- 
logical ase. 

In the present algorithm, we use the brightness at IKGHz (horizontal polarization) for the pri- 
mary rain information and 37 (ill/ to provide some help when tlx* rain does not fill a grid 3 cell 
uniformly and to enhance the sensitivity to low rain rates. The brightness temperatures as a function 
of rain rale calculated by the N\ illicit et al. ( l‘>77) model are given in f igure 10. The cell structures 
were arranged such that for each grid 3 cell there are four coincident grid 4 cells. The algorithm uses 
as input the ISM brightness temperature at grid 3 and for four corresponding (grid 4) 3711 brightness 
temperatures. Three different cases may be defined. 

1. All four 3711 brightness temperatures saturated (rain rate > 4 mm/hr). 

In this case the rain rate is derived from the I SI I brightness only using a table corresponding to 
Figure 10. 

2. All four 3711 brightness yield meaningful rain rates (rain rates <4 mm/hr). 

All 5 input v allies are converted to ram rales using the tabulated \ alues and the results are averaged. 

3. Some, but not all. of the 37 (ill/ brightness yield meaningful rain rates. 

The meaningful rain rates for each of the grid 4 cells are converted to I Kll brightness tempera- 
tures using tabulated values. Knowing these values and the average I KII brightness lor all four grid 4 
cells (i.e.. the grid 3 measurement ) an I Still/ brightness temperature and thus rain rate can be calcu- 
lated for the total area of the saturated 37 (ill/ measurements ai d the net rain rate for the whole 
grid 3 cell is calculated by appropriate averaging. A l ortran IV l’rogram which executes this calcu- 
lation is included in Appendix I. 
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SECTION V 


SIMULATIONS TO DETERMINE TM! RAIN SENSITIVITY OF 
THE GEOPHYSICAL ALGORITHMS 

In order to test these algorithms for this sensitivity to rain, brightness temperatures were calcu- 
lated according to the Wilheit et al. ( 1977 ) model lor various rain rates and a few different surface 
conditions. These brightness temperatures were then averaged with brightness temperatures for no 
ruin with various weights to simulate rain covering various percentages of the field of view. In Figure 
I I . we show the retrieved sea surface temperature as a function of rain rate averaged over the field 
of view for rain covering 1 , 3 , 10, 30, and 100% of the field of view. Figures 1 2 and 13 arc similar 
representations of wind speed and water vapor retrievals. It can be seen that the retrieval errors in- 
duced by rain become comparable to the net retrieval error due to all other causes, and thus unac- 
ceptably large, in the neighborhood of 1 .0, 0.75 and 0.5 mni/hr for sea surface temperature, sea sur- 
face wind speed and atmospheric water vapor content respectively. This is somewhat better perfor- 
mance than was reported by Wilheit et al. ( 1978). The algorithms were re-formulated, leaving out 
the 37 GHz channels, specifically to mitigate this rain sensitivity. 

This sensitivity to rain may be understood by considering two situations. If there is a nominal 
rain rate, say 1 mm/hr, which covers the field of view uniformly the rain’s contribution to the bright- 
ness temperature will be quite similar spectrally to that of non-raining clouds although much larger 
in magnitude. Thus, the retrieval algorithm for all parameters other than cloud liquid water will be 
only weakly affected If. on the other hand there is a rain rate of 100 mm/hr covering 1% of the 
field of view, although the average rain rate is the same, the situation is altogether different. The 
rain column would be nearly opaque at all frequencies and would contribute similar brightness tem- 
perature increases at all frequencies when averaged over the field of view. This apparent spectrum is 
not similar to any in the data base from which the retrievals were derived and thus is not rejected 
and can affect the retrieval of any parameter. Other rain rates and spatial distributions can provide 
a wide range of apparent spectra. The strong dependence on the spatial distribution within the field 
of view suggests great difficulty in correcting for the rain. 
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rhis in not by any means a definitive calculation. We have ignored more complicated variations 
in (Ik rain rate over the Held of view and any effect the splashing of the ram drops may have on the 
surface cmissivity. A more meaningful limit can be defined by comparing the errors between surface 
truth and retrievals from actual radiometric data with apparent rain rates also derived from actual 
spacecraft data. 
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SKCTION VI 


( ONI I USIONS 

An algorithm has been derived lor the retrieval of sea surface temperature, sea surface wind 
speed and the atmospheric overburden of water in the liquid and vanor phases. A refraction correc- 
tion for the altimeter lias also been derived. The achievable accuracies appear to be quite useful. 

For wind speeds greater than 7 m/s the wind speed precision appear to be about I meter/second and 
the sea surface temperature retrieval precision appear to be about 1 .5°(\ For low wind speeds 
t< 7 m/s) the wind speed retrieval degrades to 1 .6 m/s which is not much improvement over simply 
knowing that the wind speed is less than 7 nt/s and the sea surface temperature retrieval improves to 
less than 1 ('. 

At any wind speed, the water vapor content of the atmosphere can be retrieved to a precision 
of 0.15 gm/cm- wet path length contribution to less than 1 cm and liquid water content to about 
4 mg/cm - . 
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Table I 


Surface Models 


Wind Speeds (tn/s) 

Sea Surface Temperatures (°K) 


1 . 

271 


5 

m • 

276. 

Low Wind 

4. 

280. 


1 

6. 


283. 


8. 


286. 

12 . 


289. 

17. 

High Wind 

292. 

23. 


295. 

30. 


299. 

38. 




Tab’.- II 

Model Atmospheres 


Name 

Integrated Water Vapor 

1. U S. Standard 

1 .60 - 2.87 cm 

2. Tropical 

3.88-6.70 

3. Subtropical Summer 

4.23-6.95 

4. Subtropical Winter 

2.12-3.71 

5. Mid-Latitude Summer 

2.93-5.22 

6. Mid-Latitude Winter 

0.87-1.40 

7. Subarctic Summer 

2.07-3.11 

8. Subarctic Winter (Cold) 

0.43-0.64 

9. Arctic Winter (Cold) 

0.21 -0.32 


NOTE: Relative humidity set to 1 00"' wherever cloud liquid density is greater than zero. 


*- *• 






























Table III 


Cloud Models 


Bottom (km) 

Top (km) 

Density (gm/m 3 ) 

Max. Liquid Water (mg/cm 2 ) 

1 


0 1 

10 

1 

1 

0.3 

30 

0 

8 

0.1 

80 

7 

8 

0.2 

20 

1 

3 

0.04 

8 

1 

3 

0.08 

16 

2 

4 

0.02 

4 

6 

8 

0.2 

40 

- 

- 

0 

0 


NOTE: Density set to zero where atmospheric temperature <40°(' (limit of supercooled water). 
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Figure 1. Typical Earth Viewing Geometry Showing the Effect of Reflection of a Microwave 
Radiance Expressed as a Brightness Temperature Off a Surface such as the Ocean where T; is the 
Downwelling Brightness Temperature, Tf the Upwelling Brightness Temperature and E. R and t^ 
the Emissivity. Reflectivity and Thermodynamic Temperature of the Surface Respectively 


23 




























o 5 10 IS 20 2S 30 35 40 45 50 


FREQUENCY (GHi) 

Figure 2. I missivity of a Smooth W iter Surface at 20°C 


:4 



I igurc 3. Increase in Bright mss of the Ocean Surface (ll NUll » al 19. 35(111/ Due in Wind Speed 

(Nordberg el al ( 1 97 1 )> 
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INCREASE IN BRIGHTNESS 
TEMPERATURE PER UNIT 



FREQUENCY (GHz) 

I inure 4 Spectrum of Increase in Brightness Temperature Caused by Wind at the Ocean Surface 

(Webster et al. ( 1 *475 )> 
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ATTENUATION (DB/KM) 



FREQUENCY (GHz) 

Figure 5. Microwave Absorption Coefficient for Atmospheric Water Vapor as a Function 

of Frequency 
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. Data Processing Scheme for the SMMR Which was Launched on Moth Nimbus-Ci 

and Seasat-A 
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i igure N. Logical Flow Diagram of Statistical Retrieval Technique 
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